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Abstract. We present predictions for cosmological parameter constraints from combined measurements of second- and third- 
order statistics of cosmic shear. We define the generalized third-order aperture mass statistics <M.^ p > and show that it contains 
much more information about the bispectrum of the projected matter density than the skewness of the aperture mass. From 
theoretical models as well as from ACDM ray-tracing simulations, we calculate (Ml p ) and (M^ p > and their dependence on 
cosmological parameters. The covariances including shot noise and cosmic variance of Mj p , M^ p and their cross-correlation 
are calculated using ray-tracing simulations. We perform an extensive Fisher matrix analysis, and for various combinations of 
cosmological parameters, we predict 1-cr-errors corresponding to measurements from a deep 29 square degree cosmic shear 
survey. Although the parameter degeneracies can not be lifted completely, the (linear) combination of second- and third-order 
aperture mass statistics reduces the errors significantly. The strong degeneracy between Q. m and <x 8 , present for all second-order 
cosmic shear measures, is diminished substantially, whereas less improvement is found for the near-degenerate pair consisting 
of the shape parameter T and the spectral index n s . Uncertainties in the source galaxy redshift zo increase the errors of all other 
parameters. 
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1. Introduction 

In recent years, weak gravitational lensing by the large-scale 
matter distribution in the Universe has become an important 
tool for cosmology. Cosmic shear surveys have yielded con- 
straints on cosmological parameters without the need for mod- 
eling the relation between luminous and dark matter (bias). 
In particular, the power spectrum normal ization <xs has been 
obtain ed with less than 10 % uncertainty Jvan Waerbeke et alJ 
1200.4 . 

On the one hand, the observed sky area and thus the num- 
ber of faint background galaxies increased dramatically with 
the advent of wide-field imaging cameras mounted onto large 
telescopes. On the other hand, measurement errors have de- 
creased with further understanding of systematics together with 
new image analysis methods. These two advances were crucial 
in the evolution of cosmic shear towards a high-precision cos- 
mology tool. 

Cosmic shear is sensitive to inhomogeneities in the pro- 
jected matter distribution out to redshifts of order unity, de- 
pending on the depth of the survey. It probes scales where fluc- 
tuations have started to grow non-linearly due to gravitational 
instabilities. These non-linearities along with projection effects 
erase most of the primordial features such as baryon wiggles in 
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the power spectrum. Thus, cosmological parameters cannot be 
determined uniquely using second-order statistics alone; there 
exist substantial near-degeneracies, e.g. between Q m and cr%. 

Because these degeneracies manifest themselves in a dif- 
ferent way for shear statistics of different order, they can be 
lifted by combining e.g. second- and third-order statistics. 
An example is the reduced skewness of the convergence or 
projected surface mass density k, which has been shown to 
not, or only weakly, depend on cr 8 and thus to be able to 
break the near-degenerac y with Q m jBernardeau et al.lll997l 
Ivan Waerbeke et al Jl999h . 

Although the conv ergence cannot be observed directly, 
ISchneider et al.l dl998l) defined the so-called aperture mass 
statistics M ap , which is a local convolution of k with a com- 
pensated filter, and which can be measured directly from the 
ellipticities of the background galaxies. 

The first signi ficant non-zero third-orde r cosmic shear sig- 
nal was found by iBernardeau et al 1 ll2002ll . who measured an 
integral over the three-point correlation function (3PCF) of 
shear in the VIRMOS-DESCART survey. Fr om the same data , 
aperture mass skewness was detected later bv lPen et"al 1 J2003I) . 
and an upper limit for £!a was derived. A M ap skewness detec- 
tion at the 2cr-lev el was obtained from the CTIO survey by 
Ijarvis etalJj2004 who also derived handy expressions for the 
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In this paper, we demonstrate the improvement of cos- 
mological parameter determination from cosmic shear, using 
combined measurements of the second- and generalized third- 
order ape rture mass statis t ics. T his latter quantity was intro- 
duced by ISchneider et alJ d2004l) as the third-order correlator 
of M ap for three different aperture radii, (M ap (#i, 82, 63)} = 
<M ap (0 1 )M ap (0 2 )M a p(0 3 )). Unlike the skewness <M a 3 p (0, 8, 8)) = 
(M ap (8)M a p(8)M- dp (8)), which depends on only one filter scale 
8, the generalized third-order aperture mass statistics contains 
information about the convergence bispectrum in principle over 
the full Fourier-space. 

The reasons of employing M ap instead of the shear correla- 
tion functions are multiple: 

- M ap is a scalar quantity, therefore odd powers of it 
such as M ap have non-trivial expectation values. In con- 
trast, no scalar can be formed from tri-linear com- 
binati ons of the eight components of the 3PCF of 
shear dSchneider & Lombardill2003t ITakada & Jair3l2003t 
IZaldarriaga & Scoccimarrol2003l) . 

- From the aperture mass statistics, we get a measure of the 
residual syst ematics by its ability to separate the E- from 
the B-mode JCrittenden et all2002llschneider et all2002h . 
This is true for both second- and third-order. 

- The integral relations between (M ap ) and the bispectrum 
are much easier and numerically faster to evaluate than for 
the 3PCF dSchneider et al.ll2004 . The reason is that <M^ p ) 
is a local measure of the bispectrum, whereas the integral 
kernel for the 3PCF is a highly oscillating function with 
infinite support. 

This paper is organized as follows. In Sect. 13 we describe 
the theoretical models that are employed for the power and bi- 
spectrum. We give the definition of the second- and general- 
ized third-order aperture mass statistics and their relation to the 
power and bispectrum in Sect. [5] followed by a short descrip- 
tion of the ray-tracing simulations. Section|4]addresses the cal- 
culation of the covariance matrices of M|L, Ml„ and their cross- 

ap' ap 

correlation. Finally, in Sect.[5]we present our results on cosmo- 
logical parameter constraints from a Fisher matrix analysis. 



2. Models of the power- and bispectrum 

Statistical weak gravitational lensing on large scales probes the 
projected density field of the matter in the Universe, also called 
convergence k. All second-order statistics of the convergence 
can be expressed as functions of the two-point correlation func- 
tion (2PCF) of k or its Fourier transform, the power spectrum 
P K . Analogously, third-order statistics are related to the 3PCF 
of k\ its Fourier transform is the bispectrum B K . 

The 3-D dark matter power spectrum has been extensively 
modeled using numerical simulations. Halo model approaches 
as well as fitting formulae give very accurate descriptions of 
the quasi-linear a nd highly non-linear regime on intermediate 
and small scales dPeacock & Doddslll99a ISmith et al J [20031 
ICoorav & Shethll2002l). Throughout this wor k, we employ the 
fitting f ormula of IPeacock & Doddsl i 19961) . which was also 
used by IScoccimarro & Couchmanl J200ll) for their modeling 
of the bispectrum. 



On the other hand, the bispectrum of the cosmological dark 
matter distribution is less securely known. It is well established 
that the primordial density fluctuations were Gaussian (e.g. 
ISnergel et al.ll2003l) . In the limit of linear perturbations, they 
remain Gaussian - thus the power spectrum alone contains all 
information about the large-scale structure. However, gravita- 
tional clustering is a non-linear process and, in particular at 
small scales, the mass distribution evolves to become highly 
non-Gaussian. 

The bispectrum B K of the convergence is defined by the fol- 
lowing equation: 



{k(t x )k(h)m)) = (2n) 2 6 D (e l + t 2 + h) 

x [B K {t u t 2 ) + B K {t 2 ,h) + B^hJi)] , 



(1) 



where k is the Fourier transform of k and 6d is Dirac's delta 
function. 

We assume the field k to be statistically isotropic, thus 
its bispectrum only depends on the moduli €1,(2 of the wave 
vectors and their enclosed angle ip, B K (t\,t2) = b K (t\,(2,<P)- 
Because of parity symmetry, b K is an even function of if. 

In this w ork, we employ hyper-extended p erturbation the- 
ory (HEPT, IScoccimarro & Couchmanl l200lh for a ACDM 
Universe as a model for the bispectrum. The HEPT fitting for- 
mula fits the A^-body simulations with an accuracy of S 15 
percent, which is sufficient for our purpose. In HEPT, we can 
write 



b K (l\ , h, <P) = J] Ff\ix , h) cos m ( V ) ^ m) (A , £ 2 ), 



(2) 



with Ff = 10/7, Ff ] = hlh + titfu = 4 / 7 - The func " 
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r(m) 



tions b K are projections of the 3-D bispectrum of density fluc- 
tuations 5 which in the quasi-linear regime are given in terms 
of the power spectrum Pg. The projection is calculated using 
Limber's equation, and yields 

J ftc(w) 



xP (5 



fK(w)J \f K (w) 



(3) 



Here, /k(w) is the comoving angular distance (f K = id for a flat 
Universe) and w is the comoving distance. The lens efficiency 
function G is 

Vfhor 

3 ///of Q m f , Jk(w'-w) 
G(w) = - — -— dw p(w) , (4) 

2 V c I a(w) J /at(w') 

w 

where p denotes the probability distribution of the comoving 
number density of source galaxies. For the ray-tracing simu- 
lations we will assume that all source galaxies are at redshift 

In quasi-linear perturbation theory (PT), / (0) = = 
/ (2) — 1. In HEPT however, we have to insert for / (m \ m 
0, 1,2 the fittin g functions a, b and c respectively , as given in 
eqs. (10- 12) of lScoccimarro & Couchmanl J200ll) . These coef- 
ficients depend on the wave vector I measured in units of some 
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non-linear scale ^nl(w), the local spectral index of the linear 
power spectrum n(£) and weakly on the power spectrum nor- 
malization cr 8 and the linear growth factor. The HEPT fitting 
functions a, b and c par ametrize a non-linear generalizatio n of 
PT and were obtained bv lScoccimarro & CouchmanlfcOOll) us- 
ing A^-body simulations. In the large-scale limit, these functions 
approach unity to recover the PT results. For very small scales, 
a is constant, b and c vanish, so that the bispectrum (0 becomes 
independent of ip and thus the reduced bispectrum, which is ba- 
sically the ratio of b K and the square of the power spectrum, be- 
comes independent of the triangle configuration and takes the 
value of the hierarchical amplitude of stable clustering. 

For the sake of completeness, we also give the power spec- 
trum of the convergence, 



j dwG 2 (w)P 6 i^ 



fx(w) 



(5) 



separates the E- from the B-mode, it is an extremely useful 
tool to assess measurement errors and systematics. Moreover, 
(Mj p ) is a local measure of the power spectrum and therefore 
very sensitive to cosmological parameters. 

The ne xt-higher order quantity is the third mome nt or skew- 
ness of © JSchneider et aljl998tljarvis et alj2004 . However, 
a logical step is to generalize this statistics and allow for corre- 
lations on different filter scales 9\ , 9 2 and #3 JSchneider et alJ 



12004 . We denote this new quantity with (M^ p (9 u 9 2 ,9 3 )) = 
(M ap (di)M ll p(92)M. dp (83)), in contrast to the case of three equal 
filter scales, <M^ pd (6»)> = <A^ p (0, 9, 9)}. 

We expect the generalized aperture mass to carry much 
more information than the 'diagonal' case (M^ p d ). The latter 
basically samples the bispectrum for equilateral triangles only, 
whereas (M^ p ) probes t he bispectrum es s ential ly over the full l- 
space, as was shown in lSchneider et alJ d2004l) . see also Fig.^ 

Throughout this paper , we employ the filter functions given 
bv lCrittenden et al.N2002h . 



3. Second- and third-order aperture mass 

3.1. Definition 

The apert u re ma ss, introduced by iKaiser et alJ Jl994h and 
ISchneidell dl996l) . is defined as the integral over the filtered 
surface mass density k in an aperture, centered at some point 
&. Alternatively, it can be expressed in terms of the tangential 
shear y t {§') = -^R[y(#') exp(-2iy>)], where tp is the polar angle 
of the vector §' - §, such that the tangential component of the 
shear is understood with respect to the aperture center With 
a filter function Ug, the definition reads 



M ap (0,#) = JdVl/ e (|# -#'!)*(#') 
= j d 2 &'Q e (\#-&'\)y t (&'), 



(6) 



the second equality holds if Ug is a compensated filter function, 
i.e. / d$$Ug(&) = 0, and 
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f di?'#' 

■o 2 J 



Ug(V) - Ug{§). 



(7) 



It is the second equality in which makes the aperture mass 
statistics so useful, because it can be estimated by averaging 
over the (weighted) tangential ellipticities in an aperture. The 
integrals (J6j can be written as convolution, 



M ap (0, &) = {Ug * k) (&) = %(Q' g * y) {&), 
where we defined the modified filter function 



-2iarctan(e 2 /9i) 



(8) 



(9) 



The first moment of © vanishes, because of the compen- 
sated nature of the filter Ug. 

The second moment or dispersion of (|8j JSchneider et alJ 
has been measur ed with great success in numerous cos- 

"" 120051 IJarvis et alJ 
2003). Because it 



mic shear surveys (e.g. [van W aerbeke et al 
120031 iHoekstra etaP 120021: lHamana et all 



2n9 2 



2 f)2 



e"3JT; Q g (ff) 



47T9 4 



e 2« 2 . 



(10) 



The disadvantage of these functions is their infinite support. 
Although decreasing exponentially, they are significantly non- 
zero up to about three times the aperture radius 9. However, the 
usage of fil ter functions with finite support, e.g. the polynomial 
filters from lSchneider etal 1 ill 9981) . would involve much more 
cumbersome expressions for (M^ p ) as a function of the 3PCF, 
which makes the aperture mass statistics very unhandy when it 
has to be inferred from real data. 



3.2. Theoretical calculations of the aperture mass 
statistics 

The second- and third-order aperture mass statistics can be cal- 
culated as integrals over the power spectrum and the bispec- 
trum of the convergence k, respectively. For second order, we 
have 



—p K ({)u 2 m, 



(id 



where U{6€) = T[U e ](£) = (9{) 2 /2 ■ exp[-(6tf) 2 /2] is the 
Fourier transform of the filter function Ug. The generalized 
third-order aperture mass statistics can be written as 

(mI p (9 u 9 2 ,9,)) = <M ap (0 1 )M ap (0 2 )M ap (0 3 )) 

r d 2 fi r d% 



x J] fr(Wil) u(ej\e 2 \) u(e k \h + 



(12) 



where S 3 is the symmetric permutation group of (123), thus 
the summation is perf ormed over even permutations of j, k 
JSchneider et al.l2004 . 

Both integrals ( II It and (I12> are easily calculated numeri- 
cally due to the exponential cut-off of U for large t. Eq. d!2t 
can be simplified further if the bispectrum can be factorized as 
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Fig. 1. The filter functions l {m) (1161 for the generalized third-order aperture mass statistics as a function of the bispectrum (I15L 
Contours of the I " 1 for m = 0,1,2 from left to right are plotted for different values of 0,- as a function of l\ and l 2 (in units of 
inverse radians). The dashed contours indicate negative values. The right-most panels show the profile along the diagonal of I tm \ 
where the solid, dashed and dotted lines correspond to m -0,1 and 2 respectively. 



in (0. Then, terms of the form 

K (m \ei u eii) = 
in 

j dtp cos m (<p) ft [e ^ + 1\ + 2€\£ 2 cos <p j , 



(13) 



for m — 0, 1, 2 can be separated and carried out analytically, 



K (0 \t u t 2 ) = m -h^ 
K^(t u t 2 ) = mrW*® 

K {1 \t u t 2 ) = to"* 



[(t\ +$)lQ(hh)-2hhlx(hhj\., 



[2f 1 f 2 I (fif 2 ) 
(2 + f2 + f 2)l l(flf2 )] 

(2 + r x + tl)l (ht 2 ) 



(14) 



(' +' ) 



t ? + 1 i + 4 



fl?2 



+ 2*1*2 



Il(flf2) 



where I„ is the modified Bessel function of order n. We get for 



(M3 p (0 1 ,0 2 ,0 3 )) = (27rr 3 J d£, A j A( 2 



m=0 

with 

i (m \e u e 2 ,ey,(ij 2 ) = Ff\t u { 2 ) 



(15) 



(16) 



One sees in Fig. ^ that the functions 7 (m) as defined in the 
previous equation are relatively well localized which makes 
(Mg p ) a local measure of the bispectrum. Further, for equal fil- 



ter scales, only the region around the diagonal of the is 
probed, corresponding to equilateral triangles in Fourier space. 
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When different filter scales are taken into account, other parts 
further away from the diagonal of b™ can contribute to the in- 
tegral. Although in this latter case the amplitude of l {nt) is lower, 
the generalized aperture mass, probing the bispectrum for gen- 
eral triangles in Fourier space, contains much more informa- 
tion about the bispectrum and cosmology than the 'diagonal' 
one. This approach to sample the bispectrum o n a large region 
of Fo urier space is similar to a previous study jTakada & Jainl 
who have used all triangle configurations of the conver- 
gence bispectrum in order to predict tight constraints on cos- 
mological parameters from cosmic shear. In contrast to that 
work, we use moments of the aperture mass statistics (which 
are direct weak lensing observables) as real-space probes of 
the convergence power spectrum and bispectrum. 



3.3. Ray-tracing simulations 

We use 36 ACDM ray-tracing si mulations, kindly pro vided by 
T. Hamana (for more details see Menard et al 1 l2003l) in order 
to calculate the second- and third-order aperture mass statistics 
and their co variances (Sect. |4} . Each field consists of 1024 2 
data points in k and y, the pixel size is 0.2'. We assume our 
galaxies to be given on a regular grid - every pixel corresponds 
to a galaxy, thus our source galaxy density is 25 per square arc 
minute. We note here that the Poisson noise is much smaller 
than the shape noise of the ellipticities, and that apertures with 
radii smaller than one arc minute are discarded due to discrete- 
ness effects in the ray-tracing and in the underlying N-bo&y 
simulations. 

All source galaxies are located at a redshift of about unity. 
See Table^for the fiducial values of the parameters. 

Because the field k is given on a regular grid, moments of 
the aperture mass statistics © can be calculated very quickly 
using FFT, with the ensemble averages replaced by the average 
over all aperture centers However, since for discrete Fourier 
transforms, periodic boundary conditions are assumed, which 
is not the case for the ray-tracing simulations, points near the 
borders have to be excluded from the averaging. This leads to 
smaller effective area and therefore to an overestimation of the 
covariance of the M ap -statistics, which increases with the aper- 
ture radius. In order to avoid this, one could calculate (M^ p ) 
and (Ml p ) from the shear correlation functions, which takes 
into account the complete area. This approach is not chosen 
here because of the time-consuming calculation of the 3PCF. 
The correction scheme we apply to the covariance matrices is 
described in Sect. 14.11 



Table 1. Fiducial values of the cosmological parameters that 
are used for the theoretical model to match the ray-tracing sim- 
ul ations. If the shap e parameter T is interpreted as Sugiyama's 
r JSugivairialll995l) . our fiducial model corresponds to Q.\, - 
0.04 and h = 0.7. 
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Fig. 2. (Ml p ) and (Mj p ) from the 36 ACDM simulations 
(solid lines) as compared to theoretical predictions (dashed 
and dotted lines). The error bars are the rms values from 
the 36 fields. (M^ p ) is calculated from the simulations for 
aperture radii smaller th an one si xth of the field size. PD = 
IPeacock & Doddsl Jl996l) . HEPT = IScoccimarro & Couchmanl 
J200ll) . PT = (quasi-)linear perturbation theory. 



Figs.|2]and|3]show (Mj p ) and <M^ p ) from the ACDM sim- 
ulations and the theoretical predictions based on HEPT. The 
non-linear fitting formulae reproduce reasonably well the re- 
sults from the simulations for angular scales above ~ 1 arc 
minute. The largest aperture which can be put onto the field 
without being too close to the border is for # max = a/6 = 34', 
where a = 204.8' is the field size. 

For comparison, we calculate (ML ) by integrating o yer th e 
3PCF, using eqs. (62) and (71) from ISchneider et all J2004 
Althou gh we use the fast tree-code algorithm of I.Tarvis et all 
J2004 to calculate the 3PCF, it is still very time-taking since 
a fine binning of the 3PCF is needed (see below). Our results 
are shown in Fig. 0] and represent the average over three of 



the ray-tracing fields. (ML) as calculated via apertures can- 
not be determined for large radii because of the border effects, 
as mentioned above. Since (M^ p ) obtained via integrating over 
the 3PCF is based on the simulated shear field, we use the y 
fields instead of the k fields in order to calculate (M\ p ) via the 
FFT aperture method, using the second equality in (j8jl. With 
Mj_(#) = 3 (Q' e * yj, we also determine the statistics <M^ p M ± ), 
(M ap Mj_) and (M\) as indicators of a B-mode. (M. dp M\) is ex- 
pected to vanish if the ray-tracing simulations are B-mode-free. 
The two quantities with odd power in M ± can only b e non-zero 
for a c onvergence field which is not parity-invariant JSchneiderl 
120031) . We found all three statistics to be three and more orders 
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Fig. 4. (M^ p d )(6) from apertures (bold line), and from the in- 
tegration over the 3PCF (thin line). Also plotted are the B- 
mode signals from the integration method (M ap M^_) (dashed), 
(Ml p M ± ) (dotted) and {M\) (dash-dotted). The curves repre- 



5. 10.0 20. sent the mean from three of the ray-tracing fields. 



Fig. 3. Contours of (Ml p (Oi,02,0?,)) from simulations (upper 
row) and from the HEPT model (lower row). In each panel #3 
is fixed to the value indicated by the cross. These are 3.87 arc 
minutes (left column) and 10.77 arc minutes (right column). 



of magnitude below the pure E-mode, confirming that the ray- 
tracing simulations contain virtually no B-mode and are parity- 
invariant. However, when inferred from the 3PCF, {M. dV M\) is 
at a couple of percent of the E-mode. This is most probable due 
to the binning of the 3PCF — the B-mode gets smaller when 
we refine the binning. In our calculations, we use a logarithmic 
bin width of b — 0.075. As can be seen in Fig.[4l there is good 
agreement between the two methods, except for very small an- 
gular scales (where the B-mode is of the order 10%) and large 
aperture radii (where a significant fraction of the field near the 
border can not be taken into account with the aperture method). 

Integrating over the 3PCF is the preferred method in the 
case of real data, since the determination of correlation func- 
tions is not affected by unusable regions which makes placing 
apertures onto the observed field very ineffective. However, the 
calculation of the 3PCF is very time-consuming even using the 
fast tree-code algorithm. Moreover, a relatively fine binning of 
the 3PCF is needed in order not to introduce a B-mode from 
the in tegration of the 3P CF, and the computation time goes as 
b~ 3 3 jjarvis et alj |2004) where b is the logarithmic bin width. 
With b = 0.075, the integration method takes about a factor 
500 longer than the aperture method using FFT 




8 [arcmin] 




Fig. 5. Logarithmic derivatives of (M% p ) and (M^ d ) with re- 
spect to some cosmological parameters as indicated in the fig- 
ure legend. 



3.4. Dependence on cosmological parameters 

The goal of this paper is to study the ability of weak lens- 
ing measurements of the aperture mass statistics to constrain 
cosmological parameters. It is therefore instructive to show the 
dependence of the aperture mass on various cosmological pa- 



rameters, and to compare its second- and third-order moments. 
The more different the dependencies are for the second- and 
third-order statistics, the better will be the improvement on the 
parameter constraints when combining both. 

In Figs. 13-0 the logarithmic derivatives of the aperture 
mass statistics with respect to cosmological parameters used 
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Fig.6. Contours of dlog(M 3 dp {6i,0 2 ,0 3 ))/dp from the HEPT pig.7. Contours of dlog<M a 3 p (6» I ,6» 2 ,6» 3 ))/d/7 from the HEPT 

model, with p = Q m , Q A , T from top to bottom. In each panel model, with p = cr 8 , n s , Z q from top to bottom. In each panel 

6» 3 is fixed to the value indicated by the cross. These are 3.87 # 3 [ s nxe d to the value indicated by the cross. These are 3.87 

(left column) and 10.77 arcminutes (right column). (l e ft column) and 10.77 arcminutes (right column). 



here are shown. In all cases, the curves shown in Fig. [5] are 
quite featureless, their similarity is due to the near-degeneracies 
between the parameters. For example, we find that the ratios 
(d(Ml p )/dn s )/(d(Ml p )/dT) * (d{Ml p i )ldn s )l(d{Ml p A )ldY) are 
roughly equal and constant as a function of the aperture radius 
0. Therefore, we expect these two parameters to have the same 
near-degeneracy for both statistics. 

The ratio of derivatives with respect to D. m and erg are 
slowly increasing functions of 0, with significant differences 
between (M ap ) and (M 3 p d ). From that we can infer that the re- 
duced skewness 53 = (M 3 pd )/(Ml p ) 2 breaks the Q m - <x g degen- 
eracy of second-order cosmic shear statistics. From Fig.[3]one 
sees that dlog<M 3 pd )/<9cr 8 * 2<91og<Mj p >/<9cr 8 , so dsildo-% a 
— 53 is indeed nearly independe nt of cr 8 , as predicted 
from quasi-linear per turbation theory jBernardeau et aljll997t 
ISchneideretalJl998l) . 



4. Covariance matrices of aperture mass statistics 

4.1. Definition and Notation 

Let Mj be an estimator of some statistics, e.g. of the second- 
order aperture mass (M^ p (6j)) for some aperture radius 0j. The 
covariance matrix of this estimator is defined as 

Cov(M), 7 = (MiMj) - (Mi) (Mj) . (17) 

In case of the generalized third-order aperture mass 
(Ml p (0j, 8j, 6k)), the covariance depends on six scalar quanti- 
ties, namely the 2x3 filter scales involved. In order to ob- 
tain a two-dimensional matrix, we relabel all non-degenerate 
combinations of filter triplets {{M\ p ) is invariant under per- 
mutations of its arguments) with a single index. The result- 
ing (Ml p (6 u 6 2 , 6> 3 ))-vector is organized such that (61,62,0?,) 
is in lexical order, we further demand that 6\ < 82 < 63. 
Note that the labeling order does not play a role in the later 
analysis. For a number of distinct filter scales, there are 
^ 3 +2 j = 7V(/V + l)(N + 2)16 different combinations. 

We define the two covariance matrices Cov(M ap ) and 
Cov(M ap ) for the second- and generalized third-order aper- 
ture mass statistics, respectively. Further, for the skewness of 
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M ap , which is a function of only one filter scale, d (ff) = 
Atl p (0, 0, 0), we define the covariance matrix Cov(M^ p d ). 

The averaging in eq. dl7l is performed over the different 
simulations. Because of the small number of realizations, we 
split up each of the 36 fields into 4 subfields, corresponding to a 
survey of area A = 102.4' 2 , and average over the resulting 144 
subfields. Adjacent subfields do not represent fully indepen- 
dent realizations of the convergence field, but the correlations 
are negligible: when averaging over only a bootstrapped subset 
of subfields, we get no systematic deviation but only a noisier 
estimate of the covariance. Note that because of the splitting, 
the maximum usable aperture radius is now 17 arc minutes. 

We take into account apertures with centers not closer to 
the border than three times the aperture radius 9. This results 
in an effective area A e ff(#) which is smaller than the original 
area A - a 1 , namely A e tf(6) — (a - 6ff) 2 . Since the covariance 
is anti-proportional to the observed area, we can easily apply 
a correction scheme, and multiply each covariance matrix en- 
try Cov(6>i, fa) by y/A ef[ (8i)A eS (82)/A in the case of (M 2 ) and 
(M^L d ). For the generalized third-order aperture mass, where 
each matrix element corresponds to two triplets of aperture 
radii, the effective area corresponding to the maximum radius 
of each triplet is inserted into the correction factor. This correc- 
tion makes sure that the covariance matrix corresponds to the 
same survey area A for all aperture radii. 

For the Fisher matrix analysis of constraints on cosmolog- 
ical parameters (Sect. |5}, we scale the covariances, obtained 
from the 2.9 square degree fields, to a corresponding survey 
area of 29 square degree, by dividing them by 10, making use 
of their 1/A-dependence. Note that this increase of survey area 
is not equivalent of extending a single patch on the sky, since 
this additional observed area will not sample independent but 
correlated parts of the large-scale structure and the decrease in 
cosmic variance will be less than the increase in area. Our scal- 
ing of the area corresponds to observing 10 independent lines 
of sight, each one 2.9 square in area. 

4.2. Adding intrinsic ellipticities 

In order to realistically model the noise coming from the intrin- 
sic ellipticities of the source galaxies, one would have to add a 
random ellipticity to each shear value. It has been shown that 
this is equivalent to a dding a noise term to the convergence k 
Jvan Waerbekel2000l) . For mass reconstructions, this noise has 
to be added to a smoothed k map, however, in our case, no 
smoothing is required, thus, to each pixel of k we add a random 
Gaussian variable with dispersion cr E = 0.3/ y2. In the case 
of (Ml ), this yields the predicted amplitude to the variance 
without need of smoothing, as can be seen in Fig. [H] The shot- 
noise contribution to the v ariance is in good agreement wi th the 
Monte-Carlo method from lKilbinger & Schneider! (120041) . 

The shot-noise term of the variance of M^ p agrees very well 
with the analytical expectation JA.8L except for large 0, where 
only few apertures can be placed onto the field which are not 
too close to the border. Apparently, adding intrinsic random 
ellipticities to each grid point without smoothing introduces no 
artefacts. 



T 
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Fig. 8. The variance of (M 2 p ) from the Gaussianized ACDM 
simulations ( bold lines), in comparison wit h the Monte-Carlo 
method from lKilbinger & Schneider! d2004h (thin lines), for a 
survey area of A = 102.4' 2 . The dotted lines give the vari- 
ance from shot noise only (due to the intrinsic ellipticity dis- 
persion), the dashed curves correspond to cosmic variance only. 
The solid line include both error sources; note that it is not the 
sum of the other two curves - there is a non-vanishing mixed 
term. The dash-dotted line indicates the cosmic variance term 
for the non-Gaussian case. 

4.3. Gaussianized fields 

For Gaussian random fields, ISchneider et al 1 ll2002l) found an- 
alytic expressions for the covariance of M 2 , which were in- 
tegrated via a Monte-Carlo method by Kilbinger & Schneider 
In order to compare the results presented in this work 
with the Monte-Carlo approach as a sanity check, we trans- 
form the ray-tracing simulations into Gaussian fields without 
changing the power spectrum. This is achieved by multiply- 
ing the Fourier transform k of each convergence field by ran- 
dom phases (destroying the phase correlations). Then for each 
Fourier mode k, we pick a k(k) value randomly from one of the 
36 fields. 

Destroying the phase correlations for each individual 
field independently would not have led to the desired goal. 
Randomizing the phases cancels the connected 4-point term 
(kurtosis) of each individual realization, but not the kurtosis 
of the underlying ensemble. Our estimator of the covariance 
is independent of the kurtosis of each individual realization, 
because we first determine (M 2 p ) for each field and then av- 
erage the square of this quantity over all fields - thus for this 
averaging, only second-order quantities are taken into account. 
The process of remixing the ^-fields in Fourier space annihi- 
lates the kurtosis of the underlying ensemble, and the resulting 
fields represent realizations of a Gaussian random field. 

In Fig. [8] the variance (diagonal of the covariance) of M 2 p is 
plotted. The results from this work are in fairly good agreement 
with the iKilbinger & Schneider! d2004 Monte-Carlo method, 
although the cosmic variance term from the ray-tracings is 
slightly higher than the one from the Monte-Carlo method. 

It is clear from this figure that non-Gaussianity increases 
the noise level on the diagonal by an enormous amount, about 
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Fig. 9. Contour plots of the cosmic-variance-only term of the 
covariance of M 2 p , for the Monte-Carlo method (upper left 
panel), the Gaussianized ray-tracing fields (upper right) and the 
original fields (lower panel), for a survey area of A — 102.4' 2 . 
The contours are logarithmically spaced. 



two orders of magnitude at ~ 1'. The ratio of the non-Gaussian 
to the Gaussian variance is oc 0r 2 for small and gets less steep 
for larger 6. 

On nearly all scales, cosmic variance dominates the shot 
noise. From Fig. |9] we see that due to mode-coupling, high 
cross-correlations between different angular scales are intro- 
duced, present on the off-diagonal of the covariance. 

4.4. The case of(Ml p ) 

As for the second-order case, the variance of the aperture mass 
skewness, var(M 3 p d (6>)) = var(M 3 p (0, 6, 0)) is dominated by 
cosmic variance which is larger than the shot noise on all but 
very small scales. 

The covariance matrix of M 3 p is not diagonal-dominant, 
and shows a self-similar pattern with many secondary diago- 
nals, originating from the reordering of (M 3 p (0i, 02, #3)) into a 
single vector, which inevitably creates repeating entries of sim- 
ilar combinations of aperture radii. The correlation of (M% p ) for 
two aperture radii 0\ > 02 is a quickly decreasing function of 
the ratio 61/62- In the case of (M 3 p ) however, there are many 
combinations of filter scales which show a high correlation. 
This fact together with the small sample of realizations of k- 
fields causes the covariance matrix to be very ill-conditioned. 
For our Fisher matrix analysis (Sect. l5"2l . we have to invert the 
covariance matrix. We find stable results for the matrix invert- 
ing when the ratio of adjacent aperture radii is chosen not to be 
too small, i.e. larger than about 1.5. 



One way to determine whether our estimate of the co- 
variance of M 3 p is reasonable would involve 6-point statistics, 
which is not feasible analytically. Instead, we slightly modify 
the aperture radii used in the analysis and get a rough estimate 
of the accuracy of this method. We comment on the stability of 
our results in Sect. 15.51 

5. Constraints on cosmological parameters 

From the simulated data, we "observe" a data vector M, which 
in our case consists of the values of (M 2 ) and/or (ML) as 
a function of angular scales. Using a theoretical model, and 
approximating our observables as Gaussian variables, we con- 
struct a likelihood function £(M; p), which depends on a num- 
ber of model parameters p = (pi,p2, ■ ■ ■ , Pn)- The likelihood is 
£ oc exp(-^ 2 /2) with 

X \M-p) = J] (M k (P) - M k (po)) (Cov- 1 ),, 

kl 

x(M,(p)-M,(p )). (18) 
where the indices k and / run over the individual data points. 

5.1. The input data 

We distinguish the following five cases for the input data vector 
M and its covariance Cov: 

1. ('2') M, = (M 2 dp {0,)), Cov = covariance of M 2 p . 

2. ('3') Mi = <M 3 pd (0,)>, Cov = covariance of M 3 pd . 

3. ('3d') Mi = {Ml p {0i, 6j, 6 k )) for a combination of three fil- 
ter radii which after relabeling corresponds to index I as 
described in Sect. 14. II Cov = covariance of M 3 p . 

4. ('2+3d') Mi = some element from the concatenated data 
vector containing {M 2 p ) and (M 3 pd ). Cov is a block matrix 
containing the covariances of M 2 and M 3 , on the diagonal 
and the cross-correlation on the off-diagonal. 

5. ('2+3') Mi = some element from the concatenated data 
vector containing (M 2 p ) and (M 3 p ). Cov is a block matrix 
containing the covariances of M 2 p and M 3 p on the diagonal 
and the cross-correlation on the off-diagonal. 

Our choice of the survey geometry corresponds to ten un- 
corrected fields, each of size A = 102.4' 2 . In order to get the 
covariances of the aperture mass statistics, we split each of the 
36 ray-tracing simulations into four subfields and calculate the 
rms over the 144 resulting fields. We devide the resulting co- 
variance matrices by a factor of 10, which then correspond to 
a total survey area of A = 10 ■ 102.4' 2 = 29 square degree. We 
use six different filter radii, logarithmically spaced between 1 
and 15 arc minutes and thus have six data points for each of 
<M 2 p > and <M 3 p d ) and 56 for (M 3 p ). 

5.2. Fisher matrix 

It would be desirable to calculate the full n-dimensional like- 
lihood function in order to make predictions about error bars 
and directions of degeneracies between parameters. This, how- 
ever, is extremely time-consuming even for sparse sampling in 
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parameter space, because for every p, the bispectrum and the 
aperture mass statistics have to be calculated involving three- 
dimensional integrals. 

Instead, we use the Fisher infor mation matrix 
jKendall & Stuartlll969t iTegmarket ai1ll997l) which gives us 
a local description of the likelihood £ at its maximum. The 
Fisher matrix is defined as 

p = l d 2 [-ln£] \ = ( d 2 [-ln£] \ 
" \ dpidpj J \ dpidpj J p=po ' 

where p denotes the "true" parameter values, in our case the 
input parameters of the simulations, see Table [2 The second 
equality in dl9> holds if the maximum likelihood estimator of 
po is unbiased. The Fisher matrix is the expectation value of the 
Hessian matrix of (- In £) at p = po, which in the case of an 
unbiased maximum likelihood estimator coincides on average 
with £'s maximum - thus it is a measure of how fast £ falls off 
from the maximum. 

The smallest possible variance <x of any unbiased estimator 
of some parameter pi is given by the Cramer-Rao inequality 

<r(pd > ^J(F- l h ; (20) 

the expression on the right-hand side is called the minimum 
variance bound (MVB). 

Under the assumption that the parameter dependence of the 
covariance can be neglected, we get from eq. d 1 iS b and d!9l >: 




The derivatives of the aperture mass statistics with respect to 
the parameters pi are calculated numerically from the HEPT 
model, using po lynomial extrapolation of finite differences 
JPress et alll992l The Fisher matrix for the four combinations 
of second- and third-order statistics considered in this work is 
given in Table l5.2l 

5.3. Minimum Variance Bounds (MVBs) 

For various combinations of cosmological parameters, we com- 
pute the MVBs from the Fisher information matrix fl21i . As the 
covariance scales with A~ l (where A is the observed area), the 
MVB is roughly proportional to 1 / VA. First, the analysis is 
done for only two parameters, in order to graphically display 
the MVBs. Then, simultaneous MVBs for three and more pa- 
rameters are calculated. 

5.3.1 . Two parameters 

In Fig.^3 we show the MVBs as ellipses in two-dimensional 
subspaces of the parameter space. The hidden parameters are 
fixed. In all cases, the combination of (ML,) and (M^ p ) leads to 
a substantial reduction in the 1-cr-error. As expected, the gen- 
eralized third-order aperture mass statistics yields much better 
constraints than the 'diagonal' version (M* d ). The direction 
of degeneracy is slightly different for some parameter pairs, 
most notably when the source redshift parameter zo is involved, 



making the combination of the statistics very effective in these 
cases. The O m -crg-degeneracy is lifted partially and the com- 
bined Fisher matrix analysis yields a large improvement on the 
error of the two parameters. Contrary to that, the pair (r, « s ) is 
degenerate to a high level for both (M 2 p ) and (M;L) as well as 
for their combination. 

Note that the combined 1-cr-errors are not completely de- 
termined by the product of the likelihoods of (M 2 p ) and (Ml). 
The combined covariance is not the direct product of the co- 
variances of {M 2 V ) and (M^ p ) because of the contribution from 
the cross-correlation between both statistics. 

It is not surprising that the directions of degeneracy be- 
tween parameters are more or less similar for {M 2 p ) and (ML), 
with larger differences existing when zo is one of the free 
parameters. Both statistics depend on the convergence power 
spectrum, because in HEPT as well as in quasi-linear PT, the 
bispectrum of the matter fluctuations is given in terms of the 
power spectrum (|3}. The differences between {M 2 p ) and (ML,) 
mainly come from their different dependence on the projec- 
tion prefactor and the lens efficiency G (eq.|4}. The projection 
is most sensitive to the source redshift, and of all parameters, 
changes in zo show up in a most distinct way for {M 2 p ) and 
(M\ p ). 

Since the degeneracy directions between {M 2 p } and the 
skewness (M* d ) are very similar, not much improvement is 
obtained when these two statistics are combined and therefore, 
the corresponding error ellipses are not drawn in Fig. 1101 

5.3.2. Three and more parameters 

We calculate the MVBs for three and more parameters simul- 
taneously for various combinations of parameters and for each 
input data as described in Sect. 15.11 The results are given in 
Table [3] All hidden parameters are fixed to their fiducial val- 
ues, see Tabled If not both Q m and Qa vary, a flat Universe is 
assumed. 

In most of the cases, the error bars from the general- 
ized third-order aperture mass statistics (Ml p ) are smaller than 
those from its second-order counterpart (M 2 p ). This trend gets 
stronger the more free cosmological parameters are involved, 
since the measurement of (M^ p ) provides more data points and 
therefore more degrees of freedom 1 . The skewness of the aper- 
ture mass (M^ d ) yields by far the worst constraints on the pa- 
rameters. 

In all of the cases, the combination of {M 2 p ) and (M^ p ) re- 
sults in an improvement on the parameter constraints. This im- 
provement can be rather small, e.g. in the cases when both T 
and n s are involved. Then the combined MVB is dominated 
mainly by the MVB of (M^ p ), and the additional information 
from {M 2 p ) is unimportant. However, for a number of param- 
eter combinations, the combined error represents an improve- 
ment of a factor two and more, indicating that the dependence 
of the two statistics on the cosmological parameters is different 
to some degree, and their combination lifts the degeneracy sub- 

1 This is true only to some extent since the data points are corre- 
lated. 
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Table 2. Fisher matrix for the five different input data as listed in Sect. 15. II denoted by '2', '3', '3d', '2+3d' and '2+3', respec- 
tively. The survey is 29 square degree, all entries are given in units of 10 4 . 
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stantially. Amongst other, this occurs for the pair Q m and erg. 
Even if a rather good constraint on these two parameters from 
(Ml ) is combined with a large MVB, the combined error can 
be reduced by a factor of two and more, thus the most promi- 
nent parameter degeneracy for second-order cosmic shear be- 
tween Q. m and <x 8 can partially be broken by adding third-order 
statistics. 

When (M^p) is combined with the generalized aperture- 
mass statistics (the case '2+3') and the skewness ('2+3d'), 
the first combination always yields better parameter constraints 
than the latter. For three free parameters, the first combination 
is typically a factor of two better, if more parameters are in- 
volved, the improvement factor is even larger, up to a factor 
of ten when all six parameters are free. Thus, the preference 
of (M^p) over the skewness of M ap is justified also when it is 
combined with the second-order aperture mass statistics. 

In general, constraints on the cosmological constant £2a are 
weaker than for the other parameters, and although the combi- 
nation of second- and third order aperture mass statistics gives 
some improvement on the error, remains the least known 
parameter. 



5.4. Correlation between parameters 

The correlation coefficient of the inverse Fisher matrix 



\ " jj 

is a measure of the correlation between the / th and f h parameter. 
For ; + j, it can vary between -1 and 1. In the two-dimensional 
case, ri2 - =0 corresponds to an error ellipse with major 
and minor axes parallel to the coordinate (parameter) axes - the 
probability distribution of the parameters factorizes. For ryi — > 
1, the ellipse degenerates to a line. 

Table @] shows the correlation coefficient between all cos- 
mological parameters considered in this work. For the combi- 
nation of (M a 2 p ) and (M\ d ) ('2+3d'), the correlation is very 
large for all parameter pairs, the difference to unity in some 
cases is only of the order of 10~ 3 . The degeneracy directions of 
(M^p) and (M 3 d ) are very similar, thus the combination of the 
two causes the correlation between parameters to be very high. 
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Fig. 10. 1-cr-eiTor ellipses from the Fisher matrix. The hidden parameters are kept fixed. Dashed line: (M%„), dotted line: (M 3 ), 
solid line: (M 3 d ), dash-dotted line: combination of (M% p ) and (M^) as described in Sect. 15. II If one of the parameters is £l m , a 
flat Universe is assumed (except for the Q m -£lA-plot). The assumed survey area is 29 square degree. 



5.5. Stability 

In order to check our Fisher matrix analysis for consistency and 
stability towards small changes of the input data, we redo our 
calculations with slightly different aperture radii. For changes 
of a couple of percent in the aperture radii, the resulting Fisher 
matrix elements vary of the order of up to 10 percent. The 
MVBs (see Sect. 15. 21 fluctuate by about the same amount if 
two or three parameters are considered to be determined from 
the data simultaneously. However, for four and five free param- 
eters, the MVBs are less stable, since the Fisher matrix is nu- 



merically very ill-conditioned and the inversion is a non-linear 
operation. In general, the MVBs for {M\ p ) are less stable than 
the ones for (M% p ). 

The eigenvectors of Ft) are less affected by a different sam- 
pling of aperture radius. Angles between original and modi- 
fied eigenvectors are typically only a few degree. The variation 
of the correlation coefficient ry J22t is less than ~ 0.1 if up 
to four parameters are considered. For a higher-dimensional 
Fisher matrix however, the variation can be higher, similar to 
the case of the MVB. 
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Table 3. MVBs for various combinations of three and more 
cosmological parameters, corresponding to a 29 square degree 
survey. The hidden parameters are kept fixed. '2', '3', '3d', 
'2+3d' and '2+3' stand for the five different input data as de- 
scribed in Sect. 15. II If Q.^ is n °t a free parameter, a flat Universe 
is assumed. 
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6. Summary and Conclusion 

The power spectrum of large-scale (dark-)matter fluctuations 
was until recently the most important quantity that has been 
measured - directly or indirectly - by cosmic shear. Interesting 
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Table 4. The correlation coefficient ry ( 1221 of the inverse Fisher 
matrix '2', '3', '3d', '2+3d' and '2+3' stand for the five 
different input data as described in Sect. 15. II Note that the cor- 
relation matrix r is symmetric and unity on the diagonal. 
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constraints on cosmological parameters like Q. m and <r g have 
been obtained from second-order cosmic shear statistics. 

The bispectrum of density fluctuations contains comple- 
mentary information about structure evolution and cosmology. 
It is a measure of the non-Gaussianity of the large-scale struc- 
ture. Current cosmic shear surveys are at the detection limit of 
measuring a non-Gaussian signal significantly, and future ob- 
servations will certainly determine the bispectrum with high 
accuracy. 

Combined measurements of the power and the bispectrum 
yield additional constraints on cosmological parameters and 
partially lift degeneracies between them. The second- and gen- 
eralized third-order aperture mass statistics are local measures 
of the power and bispectrum, respectively. In this work, we 
made predictions about cosmological parameter estimations 
from combined measurements of these two weak lensing statis- 
tics. Using ACDM ray-tracing simulations, we calculated the 
covariance of M~ p and M^ p and their cross-correlation. We per- 
formed an extensive Fisher matrix analysis and obtained min- 
imum variance bounds (MBVs) for a variety of combinations 
of cosmological parameters. 

The generalized third-order aperture mass statistics 
JSchneider et al.lE()04l) is the correlator of M ap for three dif- 
ferent aperture radii. In contrast to the skewness of M ap which 
probes the bispectrum for equilateral triangles only, the gen- 
eralized third-order aperture mass is in principle sensitive to 
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the bispectrum on the complete ^-space. Therefore, it contains 
much more information about cosmology than the skewness 
alone. 

The direction of degeneracy between the cosmological pa- 
rameters considered here are similar for second- and third- 
order statistics. However, in most cases the combination of 
(M% p ) and (M% ) gives substantial improvement on the pre- 
dicted parameter constraints. The MVBs decrease by a factor 
of two or more for most of the parameter combinations. When 
the source redshift zo is not fixed but also to be determined 
from the data, the errors on the other parameters increase and 
the improvement by combining (ML) and (M^L) is lowered. 
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Appendix A: Covariance of A/.f p for shot-noise only 

Analogous to ISchneider et alJ J2002h . we analytically calcu- 



We combined the second-order aperture mass statistics late the variance of K in the case of shot-noise only, by 

integrating over the covariance of the shear 3PCF. An un- 
biased estimator of the natur al component F (0 ' of the 3PCF 
JSchneider & Lombardil2003l) is 



(M^ p ) with both the skewness and the generalized third-order 
aperture mass. The latter combination gives much better pa- 
rameter constraints than the first one. For six parameters to 
be determined from the data simultaneously, the corresponding 
MVBs are better by a factor of about 10 for each parameter. 

The Q m -cr 8 -degeneracy is very prominent for both the 
second- and the third-order statistics of M ap individually. 
However, by combining the two, the degeneracy is partially 
lifted - the 1-cr-errors of both parameters drop by a factor 
of two or more, depending on which other parameters are 
also considered to be determined from the data. The n s -T- 
degeneracy, however, can not be broken by combining (ML) 
and (Mi,), the determination of this pair of parameters is dom- 
inated by (Ml ). 

For the given range of 1 to 15 arc minutes for the aper- 
ture radii considered in this work the generalized third-order 
aperture mass statistics is dominant for the determination of 
most of the cosmological parameter combinations. The mea- 
surement error from (ML) is in general larger than the one from 
(Ma p ). However, in most of the cases, even a weak constraint 
from (Ma p ) alone contributes valuable information to the com- 
bination of the two statistics, and the combined error is much 
smaller than the one from the individual measurements. 

If the range of apertures is extended, would we expect 
the resulting improvement on the parameter estimation from 
the third-order aperture mass statistics to be higher than from 
second-order? For the former, the number of data points in- 
creases with the third power of the number of aperture radii, 
whereas for the latter, the increase is only linear. Thus, for an 
increase in the number of measured apertures, the constraints 
using (M^p) should improve more than those from (ML). On 
the other hand, the data points are not at all uncorrected; in 
fact, as it is shown in this work (Sect. I4.4L the correlation 
can be very strong for various combinations of aperture ra- 
dius triples. Moreover, for large scales (8 <; 30', see Fig. [2), 
the linear regime of the large-scale structure is probed, where 
non-Gaussian contributions are small, and the information con- 
tent of third-order shear statistics is diminished. We conclude 
that angular scales up to about 30 arc minutes will be a good 
choice for the measurement of the generalized third-order aper- 
ture mass statistics of cosmic shear. The combination of this 
statistics with (ML) will improve the resulting constraints on 
cosmological parameters quite substantially. 



ijk 



(A.l) 



where T x represents a triangle of points , # ; and uniquely 
given e.g. by two side lengths xi,X2 and the angle <p between 
them. Nj(T x ) is the number of triangles within the bin contain- 
ing T x , and S(ijk) an estimator of F (0) ; it is the following linear 
combination of products of ellipticities of three galaxies at po- 
sitions fthdj and &k, i.e. 



& — ( e ttt ~ e txx — e xtx — fixxt) + 
i( e ttx + Stxt + e xtt — e xxx)> 



(A.2) 



where = s^^&w&^&i) - s ll (#i)e v (# i )sA# i) for u, v, A e 
{'t','X'}, see eqs. (2) and (19) in lSchneider & Lombard! J2003I) . 
The summation in JA. II is performed over all possible triples 
of points ftfr), Ar z (ijk) is unity if the triangle given by 

(&i,$j,&k) is m me same bin as T x , and zero otherwise. 

The covariance of F (0) consists of four terms, which are pro- 
portional to cr\, o-g, a\ and <x", respectively. In the case of van- 
ishing cosmic variance, only the first term contributes; it reads 

Cov(f (0) ,f (0) ; r x ,r v ) = l - 

■ N T (T x )N T (T y ) 
x {& is \ijk)& {s) \lmn))A Tx (ijk)A Ty (lmn), (A.3) 

ijklmn 

where the superscript V indicates the intrinsic ('source') ellip- 
ticity. 

The term in angular brackets is non-zero only if the two 
triangles given by (#,-, §j, §k) and (§t, §,„,§„) are identical (un- 
der the assumption that different galaxies are intrinsically un- 
corrected), and factorizes into a sum of products of three 
two-point terms, each of the form (e f ,(j9 , i)e IJ f(jS'i))Sii6 l jf l i. With 
(e t e t ) = (sxS x ) = o"e/2 and (s t s x ) = 0, the term in angular 
brackets becomes 8 • [cr^/2] 3 = cr^. The sum reduces to a triple 
sum over Ajjjjk) which is just the number of triangles in the 
respective bin. Finally, we get 



Cov(f (0) ,f (0) ;r x ,r v ) = 



N T (T X ) 



S(T X , TV), 



(A.4) 



where 5(T X , 7\) is zero if the two triangles T x and T y are in 
different bins, and unity otherwise. 
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The covariance of Mj p is obtained by integrating over the 
covariance of the 3PCF. This can be done analytically for the 
case when all six aperture radii are equal (this corresponds to 
the varian ce of .) and in the ab sence of a B-mode. We write 
eq. (62) of lSchneider et alj J2004I) in the following way, abbre- 
viating the integral kernel with Rg, 

(mI p (6 u 02,ft)) = j dyi J dy 2 J # 



r(0) 



(A.5) 



where T cen denotes the 3PCF in the center-of-mass projec- 
tion, expressed as a fu nction of triangle sides as described in 
ISchneider etail (12004 . 

Before we proceed, we note that for any function /, 



j dy, J dy 2 J tyf(yuy2,iWTx,T y ) 
= f(xi,X2,<p)Ax\Ax2A(p, 



(A.6) 



where Axi, Ax2, Aip are given by the bin size in which the tri- 
angle T x is situated. 

For simplicity, we assume that boundary effects due to the 
finite field size can be neglected. Then the number of trian- 
gles within the bin characterized by (xi,X2,(p) is Nj(T x ) = 
N ■ {2nx2Ax2n) ■ {Ax\A<pri) = 2nAn i x\Ax\X2Ax2Aip, where 
N is the total number of galaxies and n is the galaxy density 
(n = N/A for A being the survey area). Thus, 

Solving the integral, we get 

/H, j 4 . 

rcmin 2 / Vlarcmin/ 



^9deg^ 



25 arc 



Note that the variance of M\ p dSchneider et aljEoO^ has the 
same dependence on the observed area A, but is only quadrati- 
cally invsere as a function of both n and 9. 
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